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Several recent studies have shown how to properly calculate the observed clustering of galaxies in a 
relativistic context, and uncovered corrections to the Newtonian calculation that become significant 
on scales near the horizon. Here, we retrace these calculations and show that, on scales approaching 
the horizon, the observed gala:xy power spectrum depends strongly on which gauge is assumed to 
relate the intrinsic fluctuations in galaxy density to matter perturbations through a linear bias 
relation. Starting from simple physical assumptions, we derive a gauge-invariant expression relating 
galaxy density perturbations to matter density perturbations on large scales, and show that it 
reduces to a linear bias relation in synchronous-comoving gauge, corroborating an assumption made 
in several recent papers. We evaluate the resulting observed galaxy power spectrum, and show that 
it leads to corrections similar to an effective non-Gaussian bias corresponding to a local /NL.eff ^ 0.5. 
This number can serve as a guideline as to which surveys need to take into account relativistic effects. 
We also discuss the scale-dependent bias induced by primordial non-Gaussianity in the relativistic 
context, which again is simplest in synchronous-comoving gauge. 

PACS numbers: 98.65. Dx, 98.80.Jk 
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I. INTRODUCTION 



The clustering of galaxies and other large-scale struc- 
ture (LSS) tracers on the largest scales has recently re- 
ceived great interest as a probe of inflation and its alter- 
natives. In the presence of primordial non-Gaussianity, 
biased tracers can exhibit a significant scalc-dcpcndent 
bias with respect to the matter distribution which in- 
creases strongly towards large scales [1]. This can be 
used as a sensitive probe of primordial non-Gaussianity 
[2]. Furthermore, ongoing, future, and proposed sur- 
veys such as the Baryon Oscillation Spectroscopic Survey 
[3], the Hobby Eberly Telescope Dark Energy Experi- 
ment [4], HypcrSuprime Cam, the Dark Energy Survey^, 
the Subaru Prime Focus Spectrograph, BigBOSS [5], the 
Large Synoptic Survey Telescope^, the Wide-Field In- 
frared Survey Telescope"^, Euclid^, and the Square Kilo- 
meter Array^ will probe modes that approach the co- 
moving horizon. All of this is strong motivation to go 
beyond the Newtonian picture of galaxy clustering widely 
adopted so far, and to embed this observable into a 
proper relativistic context. This is analogous to what has 
been done for the cosmic microwave background (CMB), 
and several aspects have long been worked out [6-8]. 
However, galaxy clustering involves a few additional com- 
plications: first, it is intrinsically three- rather than two- 
dimensional; second, one has to take into account selec- 
tion effects such as cuts on observed flux and redshift; 
and third, the relation between intrinsic fluctuations in 
the galaxy density and the fluctuations in the matter 
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density is non-trivial. 

In order to use the clustering of LSS tracers on scales 
approaching the horizon c/H{z), we thus need to under- 
stand the connection between the theoretical predictions 
for cosmological perturbations and the observationally 
inferred overdensities of galaxies. In particular, the per- 
turbations in the metric, matter density, velocity, etc. 
are always defined with respect to a particular choice of 
coordinates (gauge) , whereas observables should be inde- 
pendent of this gauge choice. 

Recently, Yoo et al. [9] have shown how to calculate 
the observed galaxy overdensity in a generally covariant, 
relativistic context. In this paper, we perform a sim- 
ilar derivation, generalizing their results in one impor- 
tant aspect. On sufficiently large scales so that pertur- 
bations are linear (and assuming Gaussian initial condi- 
tions), one commonly assumes a linear relation between 
galaxy overdensities and perturbations in the matter den- 
sity. As we show here however, the observed galaxy 
power spectrum depends on which gauge these overden- 
sities are referred to. For example, the linear bias as- 
sumed in [9, 10] is equivalent to a linear bias relation in 
the constant-observed-redshift gauge; on the other hand, 
[11, 12] adopted a linear bias in synchronous-comoving 
gauge. 

Clearly, this situation is not satisfactory, since the 
gauge choice should not impact any observable quantity. 
However, we can make progress using simple physical ar- 
guments. In a universe with Gaussian adiabatic pertur- 
bations, a galaxy knows about two properties of its large- 
scale environment: the average, "background" density of 
matter, and the local age (growth history) of its environ- 
ment. Thus, a general bias expansion should involve both 
density and age (or local growth factor). In this context, 
the simplest gauge choice is the synchronous-comoving 
(sc) gauge, where constant-time hypersurfaces are equiv- 
alent to constant-age hypersurfaces. Then, the bias with 



respect to age becomes irrelevant, and we recover the 



well-known linear bias relation: S, 



(sc) 



bS, 



(sc) 



Further 



advantages of synchronous-comoving gauge are that the 
density field in N-body as well as the output of commonly 
used Boltzmann codes are given in this gauge. We shall 
thus express most of our results in synchronous-comoving 
gauge. The transition to other gauges can be performed 
easily using expressions given in the appendices. Note 
that when properly transformed, the results derived in 
different gauges should agree. 

Recent papers by Challinor and Lewis [11] and Bonvin 
and Durrer [13] provide a further reason to reinvestigate 
relativistic corrections to the observed galaxy correlation, 
as they were not able to reach agreement with the expres- 
sions given in Ref. [9] . 

The outline of the paper is as follows. We begin by 
deriving the observed galaxy density in terms of pertur- 
bations in the synchronous-comoving gauge in Sec. II. 
In Sec. Ill, we discuss how galaxy biasing can be imple- 
mented in a gauge-invariant way. Sec. IV discusses the 
observed galaxy power spectrum. We conclude in Sec. V. 
In the appendices, we present useful results on the con- 
version between different gauges and metric conventions 
(App. A), more details on the derivations (App. B), and 
various analytical test cases for the expression for the ob- 
served galaxy overdensity (App. C). We also make the 
connection with the work of other recent papers [9-11, 13] 
in App. D. 



II. THE OBSERVED GALAXY DENSITY IN 
SYNCHRONOUS GAUGE 

In this section, wc compute the perturbations to the 
number density of LSS tracers in observed coordinates to 
linear order. 



A. Notation 

Throughout, unless otherwise noted, we adopt the 
synchronous-comoving gauge and assume a flat back- 
ground. Specifically, we write 

ds^ = a^ (r) {-dr^ + [ (1 + 2D) S,j + 2Eij] dx' dx^ } , 

(1) 
where t is the conformal time, D is a scalar metric pertur- 
bation while Eij is transverse and traceless, and related 
to the scalar perturbation E by 



E^J = [ d.d, - -<5.,V2 ) E. 



(2) 



can assign the galaxy a position x in three-dimensional 
cartesian coordinates via 



Xn, 



(3) 



where x = x(-S) and xi^) is the distance-redshift relation 
in the background Universe^ , and h is the unit vector in 
the direction of the observed position of the galaxy. In an 
unperturbed Universe, where both the observer and the 
galaxy are comoving with the background matter, this 
position in fact corresponds to the true position (in co- 
moving coordinates). This is because the photon geodesic 
in an unperturbed Universe can be written as 



x'^ix) = (to - X, fix), 



(4) 



where tq is the conformal time at observation, and we 
have chosen the affine parameter to be the comoving dis- 
tance (along the light ray). Thus, the space-time point 
of emission of the photon is given by x^^{x). 

One convenient feature of this paramctrization is that 
the geodesic equation with respect to x coincides with 
that of conformally transformed coordinates with met- 
ric gfj,iy = g^tu/a^] in the case of a spatially fiat FRW 
universe, it is simply a straight line as in Eq. (4). The 
corresponding affine parameter A in the physical FRW 
metric g^^, is determined through dX/dx oc a^. 

In a perturbed Universe, the true location of the galaxy 
is defined by the unique starting point of the geodesic 
which ends at the observer's location [xo = (0,0,0)], ar- 
rives out of the direction n, and corresponds to a photon 
redshift z (see Fig. 1); in other words, the photon fre- 
quency at arrival at Xq is 



J^o 



1 



(5) 



Here we assume we have some frequency standard i'q 
(e.g., spectral line) to compare the photon frequency to. 
In the following, it is useful to define projection oper- 
ators, so that for any spatial vector A* and tensor Eij, 

A|| =fi,X\ 
E\\ =hinjE'-^ , 



Al = {S)~fi'fij)X^, and 
Ej_ = {S'^ -h'n^)E,j. 



(6) 



Note that for a traceless tensor Eij, E± = — i?||. Corre- 
spondingly, we define projected derivative operators. 



dl=^ 



>d^ 



J' 



9|| = ffdi, and 

ai = (jy - ffhi)dj ^d'- h'd\\ 



(7) 



Latin letters denote spatial indices while Greek letters 
stand for space-time indices. 

In galaxy surveys, we observe the angular position of 
galaxies as well as their redshift z. Hereafter, we shall de- 
note observed (or inferred) coordinates with a tilde. We 



" Throughout, we assume that we have perfect knowledge of 
the background expansion history, and hence neglect Alcock- 
Paczynski-type distortions. These can be taken into account 
straightforwardly. 



Further, wc find 



dju^ = X ^{5]-n^nj), 



(8) 



from which wc derive a number of commutation relations. 
These include the commutators of the partial derivatives 
with h: 



[di,nj] = dtfij = X ^{Sij - fithj) 



(9) 



and 



[d\i,hi] = [n^dj,ni] = h^[dj,hi\ + [n^,n^]d.j = 0. (10) 

We may also find the commutators of the derivative op- 
erators with each other. Those involving the parallel 
derivatives are 

[di,d\\] = [di^h'dj] 



X ^{Sij -'hinj)dj +0 



,~,-i 



X-'d^. 



(11) 



and, since Eq. (10) shows that all components of h com- 
mute with d\\ , 



[dhd\\j\ = [d\\„d\\^] =0. 
The perpendicular derivative satisfies 



(12) 



[9x.,a||] = [d, - d\\,,d\\] = [a,,9||] = r'd^,. (i3) 

In these expressions, x is the norm of the position vector 
so that tV- = xVx- These relations are the analogue of 
the ChristofFel symbols of spherical polar coordinates. 

We also define the projection of the Laplacian operator 
(V2 = d,d') as 



9j^ = (9||,j(9^ = d\\d\\ and 



Lli^± 



X 



(14) 



Finally, we make use of 

da' = a||X|| + d^,x\_ + x^^o.ff (15) 

and 

d\\n' ^ fl'^±^ = 0. (16) 

B. Photon geodesies 

We can write the perturbed photon geodesic as 

a^^(x)=S^(x)+fe^(x); (17) 

here dx^^ is the perturbation to the photon path (note 
that the value of the affine parameter at emission is no 



x^{x + Sx) 



xf'ix) V^ ^"^^^ 




x'^ix) a(To - X) 

a{x°{x + 6x)) 



FIG. 1: Sketch of perturbed photon geodesies illustrating our 
notation. The observer is loeated at the bottom. The solid 
line indicates the actual photon geodesic tracing back to the 
source indicated by a star. The dashed line shows the appar- 
ent background photon geodesic tracing back to an inferred 
source position indicated by a circle. 



longer given by x; see Fig. 1). For our choice of affine pa- 
rameter X, we have in the unperturbed Universe [Eq. (4)] 

^ = {-hh), (18) 

dx 

while for the perturbed case [Eq. (17)] we define 

^^{-l + 5,y,h + Se). (19) 

dx 

In terms of our affine parameter x, the first-order 
geodesic equations for the fractional frequency pertur- 
bation Sv = d5x^ /dx, and the fractional perturbations 
to the photon momentum 5e' = d5x^ / dx are then given 
by 



dx 



5v^-{D' + E[{) 



(20) 



and 



dx dx\ ^ / J'^ 



D-' + iEiiY' 



E\n^ -Ewk 



(21) 



Note that the projection onto h and the spatial deriva- 
tive do not commute. In the following, we denote Eu = 



(i^ll)'*, i.e. in case of apparently ambiguous notation, the 
projection is taken before the derivative. Here, primes 
denote derivatives with respect to r, and d/dx — d\\ — dr- 
These equations are to be compared with Eqs. (9) and 
(12) in [9], respectively; note that Xthcrc = ^ Xhore and 

Before integrating Eqs. (20,21), we need to determine 
the correct boundary conditions at the observer's posi- 
tion X = 0. The observer's frame of reference is described 
by an orthonormal tetrad e" In terms of this basis, the 

components of the unit vector h of the observed photon 
are given by 



P 



(22) 



where p^ is the observed photon momentum and p = ptp^ . 
Using the metric Eq. (1) together with the orthonormal- 
ity condition 



nu a b 

9 e^e^ 



V 



we obtain 



e° =(-1,0,0,0) and 



(23) 



(24) 



where a subscript o indicates a quantity evaluated at the 
observer's position. Inserting this into Eq. (22), and re- 
quiring that h"" match the observed direction of the pho- 
ton, we obtain the following perturbations to the photon 
momentum at the observer's position (x = 0): 



Si^o=-0 and Sel== -Don' -El.n^ 



(25) 



We now integrate the spatial component Eq. (21), en- 
forcing Eq. (25) as boundary condition at x = 0: 



Se\x) 



(26) 
Integrating again up to x yields the displacements Sx^: 



2{D7V + E')fi')^ 






Pdx' 

Jo 


D^^ + E^- 


-|(n«- 


-E^fi 


{DiV + 


E^fi^o- 







6x\x) = / rfx <5e' 

= x{Dn'+E)V)o 

dx[-2{Dfi'+E)n^) 



i-X 
10 

+ (x- 



X) 



D^' + E, 



X 



{E\n'' -E«n' 



Integrating the time-component Eq. (20) of the geodesic 
equation yields 



5^{X) 



rdx{D' + E\(). (27) 

Jo 



This frequency shift contains the Dopplcr, Sachs- Wolfe, 
and integrated Sachs- Wolfe effects [6], as shown in 
App. B 1. Noting that dx'^ /dx ~ —l + 6v, we then obtain 
the time delay 



'5x°(x) 







X rx 

dx5y^- dx{x-x){D' + E'^). (28) 
Jo 



In order to obtain the perturbations to the source po- 
sition, we need to relate the affinc parameter at emission 
to the observed redshift. Recall that in synchronous- 
comoving gauge, the comoving observers' four-velocity is 
given by u^ = (a, 0,0,0). Then, the redshift along the 
perturbed geodesic at affine parameter x is given by 



z{x) 



(q '^Uf,dx^'/dx)x 
(a-2 Uf, dxt^/dx)o 

[-l + Mx)]/«(a;°(x)) 1 + ^^ 



-1 aix^ix))' 

In the second line, we have set Aq = 1 and defined 

Jz(x)^-Mx)= rdx{D' + E[). 
Jo 



(29) 



(30) 



For a given source observed at redshift z, Eq. (29) is an 
implicit relation for the affine parameter Xe at emission. 



l + z(Xe) = l + Z, 



(31) 



which defines the space-time location of the source 
through x^Q^^^j, = x'^(xe)- Note that since the confor- 
mal time of emission is r = x'^{xe)-, the redshift z{xe) 
that would have been observed for the same source with- 
out any perturbations along the line of sight is given by 
l/[l + z{xe)] = a(x°(xe)). Hence, Eq. (29) at x = Xe can 
also be written as 



l+~Z^{l + z){l+5z{Xe))- 



(32) 



To zeroth order (in the background), z(x) = ^(x)j and 
hence Xe = X- We can then expand Xe = X + i^Xi and 
Eq. (29) at first order yields 



1 + 5 = (1 + ~z)[l - {aH)i{5x^ - Sx) + 5z\ 



(33) 



Solving this for the perturbation to the affine parameter, 
we obtain 



5x = 5x° - -TTTz^Sz. 
H{z) 



(34) 



Finally, given Eq. (17), we can relate the observed posi- 
tion X, inferred assuming unperturbed geodesies x^, and 
the true position x through (see Fig. 1) 

dx 

X = X + Ax = X + x5h + -;-5x- (35) 

dx 

Separating into longitudinal and perpendicular parts, we 
obtain 



Ax II = 5x\x)fii + 5x — 5x\\ + 5x 



1 + ^A 



(36) 



and 

Axy=6x'-n'6x^i. (37) 

Eqs. (36)-(37) can be further simplified to obtain 
rx 1 _i_ 5 /-x 



Aa;||=-/ dx{D + E\\)--^ dx{D' + E\ 
Jo -"l^J Jo 



and 



Ax\ =x[E\fi^ ^E«n 



dx 



E)n^ -E\\h' 



(X - X)d\_ {D + i^ii 



(38) 



(39) 



where we have worked to first order in the displacements 



Ax. Furthermore, noting that -y, 
the background metric, we have 



-9 



where 



f 



o'5.9'; 



ff/ii/ IS 



(42) 



Finally, the galaxy density perturbations are usually 
measured with respect to the average density of galax- 
ies at fixed observed redshift, ng{z). We assume that 
when averaged over the whole survey, {5z) = so that 
{z) = z [Eq. (32)]. Also, in this paper, we follow common 
convention and define the galaxy density perturbations 
6g with respect to the comoving galaxy density. We thus 
have for the intrinsic comoving galaxy density 



a^{z)ng{x, z) = a^{z)ng{z) [I + Sg (x, z)] , 



(43) 



Note that the terms involving perturbations at the ob- 
server's location have dropped out of Eq. (38). This 
equation does not quite agree with Eq. (16) in [9], where 
i?ll has the opposite sign in the first term. This differ- 
ence also carries through to their Eq. (36). All these 
terms come from the metric perturbation dgijfi^W how- 
ever, hence they should always involve the combination 
D + E\\. For the numerical results reported in [9, 10], this 
difference is of no relevance as they evaluate the power 
spectrum in conformal Newtonian gauge where E = 0. 



C. Observed galaxy number density 

The observed number of galaxies contained within a 
volume V defined in terms of the observed coordinates is 
given by a (gauge-invariant) integral over a three-form 



N 



V-g{x°') ugix") £^^po-u^(x")tj|y fl 1. d^i 



1 



dx^ dx^ dx^ 
dx'- dx^ dx^ 



y ^"5^^") a(.TO)^'^'aii dx^ dx^ 



d^Sc 



-gugix") 



V ''^ ' a(^°) 



dx^ 



dx3 



d^x. 



(40) 



Here, x is given in terms of x by Eq. (35), e^jjypo- is the 
Levi-Civita tensor, and Ug is the physical number den- 
sity of galaxies as a function of "true" comoving loca- 
tions (in synchronous-comoving gauge). In the second 
line, we have adopted the synchronous-comoving gauge, 
where the oberver velocities reduce to u^ = (1/a, 0, 0, 0). 
In this case, not surprisingly, the volume element reduces 
to the purely spatial Jacobian \dx'/dx^\. Note that per- 
turbations enter Eq. (40) in three places: through the 
determinant \/—g', through the position- and rcdshift- 
dependence of the galaxy density n^, and through the 
Jacobian \dx'/dx^\. This Jacobian is 



dx' 



dxi 



dAa 



dxi 



= 1 



aAa 



dx' 



(41) 



where z again denotes the redshift that would have been 
observed in an unperturbed Universe, and 5g denotes the 
intrinsic fluctuations in the comoving galaxy density. Us- 
ing Eq. (32) and expanding to first order, we obtain 



a^(2)ng(x, z) = a^{z)ng{z) [1 + Sg{x)] 



(1 + ^) 



^ d{a^ng) 



dz 



Sz. 



(44) 



Note that the distinction between Sg{x.) and 5g{x.) is sec- 
ond order (this effect, analogous to CMB lensing, can 
however become important for rapidly varying correla- 
tion functions [14]). 

We can now expand Eq. (40). We define the observed 
galaxy density hg via 



/ a^{z)hg{x,z)d^x = N, 
Jv 



(45) 



so that 



a'^(z)ng(x, z) 



ff^n,(x,z-) 



dx' 



dxi 



(46) 



1 



1 + ^^.0«'(^)"9(X,2) 



1 



dAx' 



For the Jacobian, we use Eq. (15) to obtain 
9Ax* 



dSf 



:5iiAa;ii 



' ^ll^'^'ll 



■ Axwdifi' 
2Ax\\ 



-d^,Ax'j 



2k, 



X 



(47) 



where we have defined the coordinate convergence as 



\d^:Ax'^. 



(48) 



Note that the coordinate along n is defined through the 
observed redshift z. Hence, 9|| = {dxldz)d/dz. Since the 
derivative is applied to first-order displacements, it suf- 
fices to use the zeroth order expression d\\ = d / dx\x=x- 



Using Eq. (38) for Aa;||, we obtain the first two terms 
in the Jacobian 



survey depth and extinetion. Using the result Eq. (53), 
we obtain 



Axil 



\xiD^E,)-'-±!^S^ 



(49) 



and 



l + z 



9||Ax|| = - (Z? + B||) - H{z)6z — 
X dz 

Using Eq. (39) we find for the convergence 



k = K+ -V\E + -E,{o)n' - n'E'Ao), 
2 X ' 



where 



n = -\ rdx{x-x)l^l{D-\y^E + E' 
2 Jo X V 3 



(50) 



(51) 



(52) 



is the usual definition of the convergence in synchronous- 
comoving gauge and the derivatives of E with (o) are 
evaluated at the observer. The details of the derivation 
of k can be found in App. B 2. 

Finally, we can expand Eq. (46) to linear order in the 
perturbations: 



^4i!)=i + 3i? + ^, + M. + ^ 

ng{z) X 



9|l Axii — 2k 



II -^"9 



1- 



1 + S dH{z) 



H d2 



5z- 



2A.i;|| 

X 
l + z 



2k 



H{z) 



{D' + El 



(53) 



Here, ^x is given by Eq. (34), and 5z is given by Eq. (30). 
Further, we have defined 



_ dln(o 



dlna 



^_ d\R{a^ng) 



dz 



(54) 



Throughout this section, we have neglected the mag- 
nification bias contribution. This will be discussed in 
Sec. HE. 



D. Observed galaxy density contrast 



l + SdH(S) 2 l + z\ . 

-l[E'-E'{o)]-l r dx{(j> + E")-2K 
X X Jo 

-^0' + 29||i?'(o), (55) 

where (f) = D — V^-E/3 (see App. B3 for more details). 
In App. C, we apply Eq. (55) to analytical test cases 
where the exact result is known; these results serve as 
a cross-check of our result as well as to elucidate the 
significance of the various contributions. The terms in 
the first line contain the gauge-invariant intrinsic galaxy 
density perturbation (Sec. Ill) and the standard redshift- 
distortion contribution. The terms on the second line 
contain the change in volume entailed by the redshift 
perturbation. Finally, the last two lines contain further 
volume distortions from the metric at the source position, 
Doppler effect, time delay, and lensing convergence. 

The observer terms in Eq. (55), 2x~^E'{o) and 
2d\\E'{o), contribute only to the monopole and dipole of 
the galaxy distribution respectively. Therefore for most 
analyses that use the i > 2 multipoles of the galaxy dis- 
tribution, or the small angle approximation, they can be 
neglected. The monopole term is not even measurable 
since we do not know the true mean galaxy density. (The 
dipole of a galaxy distribution is measurable in principle; 
indeed it has been used to search for e.g. inhomogeneous 
initial conditions [15].) 

In App. D we connect this expression with the result 
for (Sobs of [9]. Besides notational differences we clarify 
there, the most critical difference is that we obtain a term 



Sg + beSz, 



(56) 



which takes into account the difference in the mean 
galaxy number density [Eq. (44)] between synchronous- 
comoving gauge and the uniform-redshift slicing, on 
which we measure the mean number density [the combi- 
nation Eq. (56) is manifestly gauge invariant, as we shall 
show in the next section]. This term is not considered 
in Yoo et al. [9] , as they relate the galaxy overdensity to 
the matter overdensity Sm through the gauge-invariant 
relation 



In galaxy surveys, we calculate the galaxy density per- 
turbations (5g(x) = hg{Si.)/ {fig) — 1 by referring to the 
average number density (hg) at fixed observed redshift. 
In this sense, we measure the galaxy density contrast in 
the uniform redshift gauge, where the constant-time hy- 
persurfacc is defined by Sz ~ 0. In the following, we 
assume that all fiuctuations in Ug are due to large-scale 
structure, and ignore any contributions from e.g. varying 



Sg = b{Srn -3(5z), 



(57) 



where (Sm(x) — Pmix)/pm — 1 is the fractional mat- 
ter overdensity defined in whichever gauge is adopted. 
Clearly, Eq. (57) is equivalent to assuming a linear bias 
relation in terms of density, Sg = bSm, in the uniform- 
redshift gauge. Wc will discuss these issues in the next 
section. 



On the other hand, by transforming Eq. (55) into 
conformal-Newtonian gauge, we are able to confirm that 
our result matches Eq. (30) in the recent paper by Challi- 
nor & Lewis [11]. This comparison is detailed in App. D 
as well. 



E. Magnification bias 

Equation (55) applies to the clustering of objects se- 
lected according to their intrinsic physical properties, 
their redshift, and their observed position on the sky. 
While one could in principle construct such a sample - 
e.g. a temperature-limited sample of X-ray clusters - 
most real samples in observational cosmology also de- 
pend on the apparent flux from the source. That is, they 
have a selection probability that depends on how the lu- 
minosity distance L'l differs from the mean luminosity 
distance Z)l(z). The sample selection may also depend 
on the angular size of the source, but this is not indepen- 
dent since the conservation of photon phase space density 
relates the angular diameter and luminosity distances 



i?L = (1 + zfDp 



(58) 



This section evaluates the additional terms that appear 
in Eq. (55) in the presence of a dependence on magnifi- 
cation. The key parameter that we need to measure is 
the magnification 



M 



Da 



d; 



2/'2\ 



D^^iS) Dl\~z) 



(59) 



which has mean value 1 and represents the perturbation 
to the solid angle or flux of a source, relative to a source 
at the same observed redshift z in the unperturbed Uni- 
verse. We may also write 5Jv[ = M. — 1. The galaxy 
overdensity is then 



(5g(no mag) + Q5M, 



(60) 



where <5g(no mag) is the overdensity computed from 
Eq. (55) and 



Q = 



d\mi„ 



d\-aM 



(61) 



is the dependence of the observed number counts on mag- 
nification. For a magnitude-limited sample with cumula- 
tive luminosity function n{> L) we have Q — ~d\i\n{> 
L)/d\nL. A more general criterion, e.g. one that in- 
cludes a "size" cut to reject stellar contamination [16], 



^ In the presence of opacity of intergalactic medium (IGM), e.g. 
due to Thomson scattering, this is not necessarily true. However, 
IGM opacity is a very small effect and we do not consider it in 
this paper. 



would have Q that must be determined by simulating 
the observations. 

From Eq. (60) we see that we need only determine how 
SM depends on the metric perturbations in order to have 
a complete description of the magnification bias. 

Fortunately, we have already constructed the key in- 
gredients in evaluating Ai . If we consider a right-handed 
3-dimensional orthonormal basis {n^ , a"^ , P^} where n* is 
the direction of observation, then the angular diameter 
distance can be inferred from the area perpendicular to 
the line of sight spanned by the rays along the past light 



cone near n 



D\ = V-gix") e 



fii^pa f- 



't 



dx'' dx" 

— :r~ :r~"( 



(62) 



where the partial derivatives are taken at fixed x and l^ is 
a unit purely spatial vector {u^i^ = 0) pointed away from 
the observer (in the sense that a photon emitted from the 
source with 4- velocity parallel to the null direction u'^—i'^ 
reaches the observer) . Using the chain rule to replace the 
derivatives with those involving i* gives 



TW"^ = 



D2 



V-g(^") 

HxW 



' fiupa 



Ox^ Ox^ 



(63) 



Next we observe that f is parallel to the spatial part of 
L^ = dx'^ /dx\in since it is the spatial part of the tangent 
vector to the past light cone. Since & is a unit vector 
and L'^ is a past-directed null vector, it follows that 



using Ua = (—a, 0, 0, 0) and u^ = we find 
nk L'' 1 dx'^/dx 



i{xO)L° a{x°)~dx°/dx 



(64) 



(65) 



Using this and collapsing the Levi-Cevita symbol to 3 
dimensions, we find 



M' 



V-9{x") 



^abc i^~ 



dx°- dx'' dx'' i,, ^, ~, 

—n a P^ . 



[a{x)a{x^)Ydx^ / dx " '^ dx'^ dx^ dx^ 



(66) 

Since {n\ a^ , 13^} form an orthonormal basis, we simplify 
this to 



M'^ = 



V^gi^ 



[dix)a{xO)]^dxydx 



dx^ 



dx3 



(67) 



Now we are in a position to c ompute the pieces of 
Eq. (67). We already know that y^-g{x") =a^{l+ 3D) 
and \dx'^/dx^\ = 1 + dAx^/dx'^ [the latter is given by 
Eq. (47)]. Finally the null condition gives 



dx° 
dx 



(1+D + E« 



dx\\ 



l + D + En+dnAxn. (68) 



Wc thus find 



M- 



a(a;0) 



a(x) 



1 + 2D - En 



Aa 



or 



Aa;|| 

5M = -26z -2D + En - 2 ^ 

X 



2k 



2k. 



(69) 



(70) 



This makes sense: the perturbation to the magnification 
contains the obvious coordinate convergence term 2k, but 
it also has three other pieces: a contribution —2lS.x\\/x 
associated with bringing the source closer to or farther 
from the observer; a contribution 2{—5z — D) associated 
with the isotropic conversion from coordinate distances 
to physical distances (itself having both a part from the 
change in scale factor at the source and the metric per- 
turbation); and a part E\\ associated with the anisotropy 
of the coordinate system [Eij is traceless and magnifi- 
cation depends only on the perturbation to transverse 
distances). 

Expanding Axji using Eq. (49) and k using Eq. (51) 
gives an alternate expression 

5M = -25z-2<f)--\dnE-dnE{o)]+2K-2dnE'{o) 

X 



;(^ + E||KY+i^ 



5z 



(71) 



The integral may be simplified by replacing D -\- E\\ -^ 
(j) + d'uE and then doing a double integration by parts 
using Eq. (B26), 



SM = -2(P+-[E' -E'(o)]+2K-2dnE'(o) 
X 



2 /■* 

-^ / {c^ + E")dx- 
xJo 



„ 2 1 

-2+- 



xH{~z) 



Sz. (72) 



III. 



GALAXY BIAS IN A RELATIVISTIC 
CONTEXT 



In order to make progress from Eq. (55), we need to 
relate the intrinsic galaxy overdensity 6g (here written 
in synchronous- comoving gauge) to the matter and met- 
ric perturbations. Fortunately, since we are interested in 
large scales, we only need to consider terms linear in per- 
turbations. What are the relevant quantities on which 
the physical galaxy density might depend? 

The most important characteristics of the large-scale 
environment of a given galaxy are its mean density, and 
the evolutionary stage (proper time since the Big Bang, 
or linear growth factor). In fact, on sufficiently large 
scales, these are the only quantities of relevance to the 
galaxy two-point correlations [17] ®. 



We can formalize this statement by considering some 
large spatial volume within the Universe centered around 
the spacetime point x!^, on a constant-age hypersurface, 
tu = constant, where tu denotes the proper time of co- 
moving observers since the Big Bang. Then, the number 
of galaxies (or, more generally, tracers) within that vol- 
ume can only depend on the enclosed mass M , and the 
age of the Universe in that volume tu which is being kept 
fixed: 

Ng = F(M; tu; x>^). (73) 

Here, the explicit dependence on x^ indicates any 
stochasticity in the relation between Ng and the local 
density and age. We now assume that the volume V is 
large enough so that linear perturbation theory applies. 
Then, in a given coordinate system (t, x), the enclosed 
mass is given by 

M= / p = p„,{t)[\ + 5\Yip\V 
Jv 

= PmiT)[l + 5,„-3aH5T]V. (74) 

Here, Pm{T) is the average (physical) matter density in 
the background (equivalent to p averaged over the entire 
constant-coordinate-time hypersurface), while Sm is the 
matter density perturbation on a constant- coordinate- 
time hypersurface. The second line follows from Ship = 
p/pm — 1 and dlnpm/dT = —3aH, and wc have defined 
— (5r(x) to be the displacement in coordinate time corre- 
sponding to a tu ~ constant hypersurface: 



a(T)[T — (5r(x)] = tu = constant. 



(75) 



Thus, the term —SaHSr in Eq. (74) comes in from 
going from a constant-age hypersurface to a constant- 
coordinatc-time hypersurface. Note that in Eq. (74) the 
perturbations are to be considered averaged over the vol- 
ume V. 

In exactly the same way, we can define the aver- 
age (physical) galaxy number density fig on constant- 
coordinatc-time hypcrsurfaces. The same reasoning lead- 
ing to Eq. (74) yields 



N„ 



Ug — ng(r)[l -t- 6hmg]V 



ng{T) [l + Sg-\- bep aHSr] V, 



(76) 



where bep = dlnng/dlna. We can now equate this to our 
general ansatz Eq. (73), 

Ng^FiM;tu;x^) 

= F{pmV; tu)[l + b{5m - 'iaH&T) + e] 

= ng{T)[l + 5g + bep aHSr] V. (77) 



This assumes that there is no orientation-dependent selection of 



galaxies. Such a selection will introduce a dependence on the 
large-scale tidal field as well [18]. 



Here, we have defined F{M, tjj) 
introduced tlie bias 



{F{M,tu,x^))t„, and 



^^ dlnF{M;tu) 



dlnM 



dlnF{pvV-tu) 



PmV 



d\np\ 



and the stochastic contribution to galaxy density 



e{x^') 



F{M,tu,x'') 
F{M,tu) 



- 1. 



(78) 



(79) 



Also pv denotes the average matter density (on the 
tjj — const slice) within the volume V. 

At first order, the bias b defined in this way is only a 
function of r. The stochastic contribution e to galaxy 
clustering is only a function of the spacetime point (e is 
here considered to be first order as well). In the back- 
ground {Sm — > 0,(5t — > 0), Eq. (77) implies, not surpris- 
ingly, ng{T)V — F{pmV;aT). To first order in the per- 
turbations, recall that Eq. (77) must hold in any coordi- 
nate system. Thus, we conclude that the galaxy density 
perturbation is given in general by 



Sgixf^) = b{T)[6^{xn - 3aH{T)dT{xn] 
-b,piT)aH{T)ST{x^)+e{xn- 



(80) 



On sub-horizon scales, aHdr becomes negligible com- 
pared to S,n (for standard choices of gauge). Eq. (80) 
shows that in this limit we recover the usual bias rela- 



tion Sn 



i+e. Furthermore, if we choose synchronous 



gauge where all comoving observers are synchronized so 
that tu = ar everywhere and thus St = 0, the linear 
bias relation holds on all scales. Note that the definition 
Eq. (78) is precisely what is commonly called a peak- 
background split bias parameter [19, 20]. 

This derivation was phrased in terms of the physical 
galaxy density. The reasoning and Eq. (80) trivially hold 
for the comoving galaxy density a'^ng as well, the only 
difference being that dg is now the fractional perturbation 
in comoving number density, and b^p is replaced with 
Eq. (54), 



6.= 



dln{a'^ng) 
dlna 



From now on, we shall exclusively consider comoving 
number densities, as we did in Sec. II. 

The bias relation Eq. (80) holds in all gauges, and the 
bias parameters b and b^, do not depend on the gauge 
choice. This is not very surprising, since these parame- 
ters are in principle observable: b [Eq. (78)] quantifies the 
response of the galaxy number in a given volume at fixed 
age of the Universe to a change in the average mass den- 
sity (or enclosed mass) within this volume; be [Eq. (54)] 
quantifies the dependence of the average (background) 
number density of galaxies on the age of the Universe. 
To see the gauge-invariance of these bias parameters ex- 
plicitly, consider the effect of a change in the time coor- 
dinate. 



T, 



(81) 



where T can in general be a function of r and x (spatial 
gauge transformations do not affect the density pertur- 
bations at linear order). By using Eq. (A3) in App. A, 
we find that S,-, 



and Sg transform as 



^m ^'n 



^9=^9- 



3aHT and 

^T ^5„-be aHT. 



din (a'^ng) 

dr ' ^'^^ 



(82) 



On the other hand, e is gauge- invariant. Note that a 
change in time coordinate (slicing) implies a change in 
the redshift perturbation Sz [Eq. (32)] through 



(H-z)(5z = {l + z)Sz + HT, 



(83) 



with dz/dr = —H. 

The previous two equations clearly show that the 
combination 5g + beSz appearing in Eq. (55) is gauge- 
invariant, independent of any bias relation. On the other 
hand, under the same gauge transformation with fixed b 
and be, Eq. (80) changes as 



Sg = biS,n - SaHSr) - be aH{5T + T) + e 
~ 5a — be aHT, 



(84) 



where we have used Eq. (54) in the second line. That is, 
we recover the gauge transformation of Sg with fixed bias 
parameters, and, in this sense, the bias parameters de- 
fined through Eq. (78) and Eq. (54) are gauge-invariant. 
We now see that a gauge-invariant expression for the 
galaxy number density [Eq. (55)], and a gauge- invariant 
definition of the galaxy bias are separate issues. In 
Refs. [9, 10, 13], the second issue was not addressed ex- 
plicitly, and in Refs. [9, 10] be was implicitly set to zero. 



A. Bias parameters from universal mass function 
approach 

In this section, we show how both b and be can be 
estimated in the universal mass function approach. We 
adopt the synchronous-comoving gauge, which is implicit 
in the reasoning of this approach. The universal mass 
function approach is expected to be valid for objects se- 
lected via a proxy for halo mass; however if the selection 
criteria are sensitive to merger history (e.g. one selects 
active galactic nuclei) then the universal mass function 
may not be valid. This is analogous to the merger bias 
effect in models with primordial non-Gaussianity [2, 21]. 

In this picture, one assumes that galaxies form inside 
density peaks in Lagrangian space whose height exceeds 
some critical matter density contrast 6c (i.e., p > (1 -I- 
5c)Pm)- In regions with large-scale overdensity 5i, this 
threshold is effectively lowered to 5c — 5i. If we denote 
the average abundance of tracers of mass M as n{M, 5c), 
the galaxy density contrast on large scales is then linearly 
related to the matter density contrast via 



5g{M;5i) 



1 



d\-Dn{M,5c,T) 

Wc 



5i = bSi, 



(85) 
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if wc truncate the Taylor expansion at linear order. The 
first term comes from mass conservation when transform- 
ing from Lagrangian to Eulerian space. As shown in the 
previous section, this argument is not in general correct 
in the context of general relativity. Note that if we were 
to choose a non-synchronous gauge (such as conformal- 
Ncwtonian or uniform-rcdshift gauge), the density field 
in different regions would be at different evolutionary 
stages, so that the collapse threshold Sc is not simply 
a constant on a constant-coordinate-time hypersurface. 
Therefore, the galaxy density contrast must also depend 
on the evolutionary stage, or age of the universe, in the 
region considered. However, in a synchronized gauge 
where St = 0, the argument leading to Eq. (85) is ap- 
plicable. Thus, Eq. (85) is a valid bias parameter which 
can be used in the correct, gauge-invariant bias expansion 
Eq. (80). 

Wc can also obtain a useful analytical estimate for be, 
assuming that the abundance of galaxies follows a uni- 
versal mass function, 



n{M) 



A'P 



/(^) 



diner 



dlnM 



First, the linear density bias is given by [Eq. (85)] 

dh\n dhif{v) 1 

dSc dv (T 

On the other hand, b^ is given by 

d\n{a^n) dlnf{iy) dlnf{iy) dv da 



(86) 



(87) 



dhia d\\\a dv da dlna 

V \ dhia 



(1-5). (-9. 



d\\\a 



(88) 



The logarithmic derivative of a can be further simplified 
via linear perturbation theory as 



diner d\i\D 



dlna d\aa 



= /• 



(89) 



Here / « ^^^ is the usual logarithmic growth rate fa- 
miliar from redshift-space distortion theory [22]. In sum- 
mary, we obtain 



be^{b-l)avf^6J{b-\). 



(90) 



Note that the abundance of rarer, more strongly biased 
halos evolves faster (larger 6e), and that the overall rate 
is set by the growth rate /. Eq. (90) is useful for esti- 
mating the magnitude of the corrections to the galaxy 
power spectrum. Note however that the universal mass 
function prescription might not be a good description of 
actual tracers whose redshift evolution is influenced by 
non-gravitational physics (such as star formation, feed- 
back, rcionization, etc) . The key point however is that for 
any given survey, be is in fact observable, if the redshift- 
depcndence of the source selection function is known. 



B. Bias in synchronous-comoving gauge 

In the synchronous-comoving gauge assumed in our 
derivation in Sec. II, there is no perturbation to the 00- 
component of the metric. Moreover, the constant-time 
hypersurfaccs are orthogonal to the velocities of comov- 
ing observers (in other words, v = Q). In this gauge, every 
comoving observers' proper time is synchronized, and all 
observers on a given r =constant hypersurface are at the 
same evolutionary stage, which implies 5t = 0. The den- 
sity field in standard A^-body simulations is also defined 
precisely in this gauge [23] . 

Wc now sec that the gauge-invariant bias relation 
Eq. (80) is equivalent to the well-known linear bias re- 
lation between 5g and 5m in synchronous gauge. 



A(sc) 

"9 






(91) 



where the superscripts denote that the variables are de- 
fined in synchronous-comoving gauge. Inserting this into 
Eq. (55) then yields the observed galaxy overdensity, 
completely described by the metric and matter perturba- 
tions and two numbers specific to the tracer population: 
the linear bias b and the count slope be- Note that the 
latter parameter is observable in galaxy surveys, while 
the bias 6 is a parameter that needs to be fitted for. 



C. Primordial non-Gaussianity 

In the picture outlined in this section, it is also 
straightforward to understand the effect of primordial 
non-Gaussianity. Consider a constant-age hypersurface 
as defined earlier in this section, at some early time long 
before the tracers of interest formed. Since the linear 
growth factor is the same everywhere on this slice, the 
variance of the small-scale density field a\ smoothed on 
some scale R is also the same everywhere, in the case 
of Gaussian initial conditions. In the presence of non- 
Gaussianity of the local type, mode-coupling induces a 
modulation of cr|j by long-wavelength (Bardcen) poten- 
tial perturbations ^i, so that, within a region on a 
t[/ ~ const hypersurface where <f>L can be considered 
constant, it is given by 



o-|f=-cr^(l + 4/NL*L). 



(92) 



Here, cr|j is the variance derived from the Gaussian part 
of (f>. We see that this is closely related to perturbing 
the local age of the Universe, which leads to a change in 
the local a\ as well. The relation between $ and 5m in 
synchronous-comoving gauge is given by (e.g., [27]) 



(5^'=)(fc,z) = 7W(fc,z)$(fc,z,) 
^ 3anggfc^(l + z) 
2T{k)D{z) 



$(fc,^,), (93) 



where z^, is some reference redshift where the non- 
Gaussian parameter /nl is defined (for example, that 
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of the last-scattering surface) . Using the universal mass 
function prescription (Sec. Ill A), we then see that the 
bias relation Eq. (80) in synchronous-comoving gauge is 
modified to 



9 



[b + 2h^{b-l)5,M-\k)\5t'\ (94) 



in agreement with [1]. It is straightforward to generalize 
this derivation to more general types of non-Gaussianity 
[27, 28]. Note in particular that (for the local case), the 

scale-dependent correction to 5g is proportional to fc~^ 
out to arbitrarily large scales (sec also [12, 29]). 



IV. THE LARGE-SCALE GALAXY POWER 
SPECTRUM 



We now calculate the observed galaxy power spectrum 
including the bias relation and the volume effect we have 
calculated in the previous sections. Throughout this sec- 
tion, we use the cosmological parameters from Table 1 
("WMAP-HBAO-hi/o ML") of Komatsu et al. [24] as our 
reference cosmology. As explained in App. B3, we ne- 
glect the lensing contribution k, as it is not simply incor- 
porated into a three-dimensional power spectrum. Fur- 
ther, we neglect two very small contributions, the inte- 
grated Sachs- Wolfe (ISW) contribution to (5z, and the 
time-delay contribution ex / dx (0 + E")/x. We also ne- 
glect the stochastic contribution e to 5g in the following, 
and assume that the primordial density perturbation fol- 
lows Gaussian statistics. 

Neglecting the ISW contribution, we have (App. B 1) 



5z{x)^d\\E'{x) + E"{x). 



(95) 



where we have dropped the unobservable, constant con- 
tribution from the perturbations evaluated at o. Using 
Eq. (55), Eq. (72), and the results from Sec. Ill, the 
observed galaxy density contrast written in terms of per- 
turbations in synchronous-comoving gauge is then given 
by (see App. B for the derivation) 



hSt 



be{d\\E' + E") + 2{\-Q)(j}- 



^lE' 
aH 



where 



-^il-Q)E' + {C-l)idiiE' + E"), (96) 



1 + zdH 1 + z2 , ^, ^ ,, 

--(1-Q)-2Q. (97) 



H dz 



H X 



In a ACDM Universe, the first term in C can be simplified 
to yield 

C = ln^{z)-^l{l~Q)-2Q. (98) 

2 H X 




FIG. 2: Coefficients of new terms in observed galaxy overden- 
sity in synchronous conioving gauge Eq. (100) introduced by 
relativistic volume effect and bias. Here, we plot for 6 — 1.5 
(solid), 6 = 2 (dotted), 6 = 3 (dashed) and 6 = 4 (dot-dashed) 
cases, and 6e is calculated assuming the universality of the 
mass function [Eq. (90)]. We ignore the magnification effect 
by setting Q = 0. 



where f7m(z) is the matter density parameter at redshift 
z. In Fourier space, Eq. (96) reads 



.PE' 
aH 



dg{k) = hb^^f^ + h,{ikiiE' + E") + 2(1 - 2)0 ■ 

-(l-Q)^ + (C-l)(^fcM^' + ^"), (99) 

X 

where /i is the cosine of the wave-vector fc with the line- 
of-sight direction. 

Finally, by relating the synchronous-comoving gauge 
metric perturbations to the density contrast bm 
(App. A3), we can further simplify bg as 



l^=&+/M^+4+^^, 



(100) 



where x = hj aH is the wavenumber in units of the co- 
moving horizon, and we have defined the coefficients 



J\ — c\ r^ 



2/ 



2f 



and 



S=/[6e+C-l] 



22 
(101) 

(102) 



For given value of 6e and Q, A and B are only functions 
of redshift, and incorporate the relativistic bias as well 
as the volume distortion and magnification effects. On 
small scales, when x ^ 1, we recover the usual Fourier- 
space galaxy overdensity in "Newtonian theory" , 



"9 



[h + ![i^)b^ 



(103) 
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FIG. 3: Same as Fig. 2, but for fixed bias (b = 2) and varying 
magnification with Q = — 1 (solid), 0.5 (dotted), 1 (dashed), 
2 (dot-dashed). For diffuse backgrounds, Q = 1, and the 
magnification effect cancels almost all volume distortions. 



In principle, be can be measured from the survey it- 
self provided one has good knowledge of the redshift- 
dependence of the selection function. In the following, 
we will use 



SJib~l) 



(104) 



as predicted by the universal mass function ansatz 
[Sec. Ill A] for our illustrations. Figure 2 shows A and B 
as function of redshift when we only include relativistic 
volume and bias effect, i.e. Q = 0. Note that both A and 
B diverge to — oo as z — >■ 0, because of the 1/x factor in 
C [Eq. (97)] which is a result of the volume distortion by 
velocities and gravitational redshifts. Figure 3 shows A 
and B for fixed bias (6 = 2), and varying magnification 
coefficient (Q = -1, 0.5, 1, and 2). For Q = 1, which 
is the case for diffuse backgrounds, magnification effect 
cancels volume distortion, and both A and B are small. 
For sufficiently high redshift when the Universe is approx- 
imately matter dominated, C -^ 317^/2 — 2Q ~ 1.5 — 2Q 
and / — > 1, so that 6e — >■ (5c(fe — 1) and we can approx- 
imate the coefficients as .4 — > 5.2 — 6Q + 6c{b — l)/2, 
B^0.5-2Q + Scib~l). 

From Eq. (100) we can calculate the observed galaxy 
power spectrum in terms of the linear matter power spec- 
trum P'^^'^^k) in synchronous comoving gauge, yielding 



P(sc)(fc) 



(b+f^'r+2{b+f^')^+-,+ 



A 

2 



A^ fi^B 



2r2 



(105) 



Note again that we have neglected all projected quanti- 
ties here, most importantly the magnification contribu- 
tion —2k. Furthermore, the flat-sky calculation employed 
here is likely not applicable to the very largest scales in 
actual galaxy surveys. Here we are mainly interested in 
the issue of galaxy biasing however, and defer the calcu- 
lation of the full angular galaxy power spectrum to future 
work. 
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FIG. 4: Three different theoretical predictions of the observed 
galaxy power spectrum on large scales for galaxies with linear 
bias 6 = 2 at redshift z = 0.5, and assuming Q = 0: Newto- 
nian linear theory [22] (dotted line), relativistic linear theory 
with linear bias in uniform redshift gauge [9] (dashed line), 
and relativistic linear theory with linear bias in synchronous 
comoving gauge (this work, solid line). We show both line- 
of-sight directional power spectrum (/i = 1, thick lines) and 
perpendicular directional power spectrum (/i = 0, thin lines). 
The vertical solid line indicates k — aH ai z — 0.5. 



Fig. 4 shows the galaxy power spectrum in three differ- 
ent calculations, each for /x = and fi = 1. The dotted 
lines show the linear small-scale limit given by the Kaiser 
formula. The solid lines show the relativistic calcula- 
tion Eq. (105) using the galaxy bias prescription derived 
in Sec. Ill, i.e. a linear bias relation in synchronous- 
comoving gauge. Clearly, the prediction departs from the 
Kaiser formula on scales k < 10~^ /i/Mpc. The dashed 
lines in Fig. 4 show the result of [9] for comparison. As 
we have seen in Sec. Ill, their result is equivalent to linear 
biasing in the uniform-redshift gauge. We see that the 
departures from the small-scale limit are much more sig- 
nificant in this latter calculation, showing that the precise 
choice of bias relation is important on very large scales. 

Fig. 5 shows the two-dimensional galaxy power spec- 
trum from Eq. (105), illustrating the evolution of the 
angular dependence of Pg{k,fi) with scale. Also shown 
for comparison is the prediction of the Kaiser formula 
[22]. Again, deviations appear for k < 10^"^ /i/Mpc and 
become more significant for transverse separations. 



A. Effective primordial non-Gaussianity 

Our results show that the observed power spectrum 
departs from the small-scale calculation on sufficiently 
large scales, with terms proportional to (k/aH)^^ and 
[k/aH)^^. This is reminiscent of the scale-dependent 
bias induced by primordial non-Gaussianity [1], see 
Sec. inc. It is thus natural to compare the two effects. 
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FIG. 5; Observed galaxy power spectrum (Eq. (105)) as function of feji and k±. The color contours and dashed lines show 
the Newtonian result from Kaiser [22] , while the solid lines show the result including relativistic corrections. Here we have set 
Q = 0. For reference, we also show the real space power spectrum contours as white dotted lines. 



Neglecting all relativistic effects, and non-Gaussian 
redshift distortion terms that are unimportant on large 
scales [25], the observed galaxy power spectrum for local 
primordial non-Gaussianity is given by 



3(NG) 



(fc) 



PM(fc) 



{b + Abik) + ftiy 

{b + ffi'-'f + 2Ab{k)(b + ff,^) 



Here we have defined the bias correction 



Ab{k) = 



2/NL(&-l)(^c3a„ogg 
2D{z)k^T{k) 



/nl4(&-1) 



D{z)k^T{k) 



A^™(z)<5c(6-1)/nl;^, 



(106) 



Ab^{k). 



(107) 



where in the last line we have assumed that k < 
0.01/i/Mpc so that T{k) = 1. We can define an ef- 
fective non-linearity parameter f^ such that the scale- 



dependent non-Gaussian bias inserted into the small- 
scale expression Eq. (106) leads to a power spectrum 
matching the actual observed power spectrum in the 
Gaussian case including the relativistic effects. How- 
ever, the relativistic corrections affect modes perpendic- 
ular and parallel to the line of sight differently, while the 
non-Gaussian scale-dependent bias is isotropic. [Note in 
particular the B term in Eq. (105).] Therefore, we can, 
in principle, distinguish the relativistic effect from the 
non-Gaussian scale dependent bias through the angular 
dependence of the 2D power spectrum. Given the limited 
number of modes on scales k < 10~'^ft./Mpc, it is more re- 
alistic to consider the effect on the angle-averaged power 
spectrum monopole 



pt'Hk) 



d^iPg{k,n). 



(108) 



From Eq. (105), we find for the relativistic corrections to 
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the monopole power spectrum 



j(e=o) 



(fc) - Pc 



(Ka 



,i=0) 



(k) 



P(sc)(/j) 



1 



A 



^ + ^, (109) 



a;' 



and from Eq. (106), we find that the non-Gaussian cor- 
rections only are given by 



j(NG,«=0) 



P/^"''="Hfc) - Pg^ 



(Kaiser, f=0) 



(k) 



P(sc)(fc) 

-^f2m(5c(6 - 1)/nl 



X 



f 



where 

p(Kaiscr,£=0)/^\ 



■Diz) 



'^^^\M^\f 



Ab^k), 



(110) 
Pt\k) (111) 



is the monopole galaxy power spectrum given by the 
Kaiser formula. By equating the term proportional to 
x~^ in Eq. (109) and Eq. (110), we find the effective am- 
plitude of non-Gaussianity as 



Dia) 



1 



rcS 

■'^^ a 6n,n6c{b-l) 



2A + 



B^ 



3b + f 



(112) 



I Dja) 
2 a 



/I 



2/ 



l + 2//a„ + C-/-2Q 
Sc{b-1) f,. 



3n„,i3b + f) 



r + fT 



C-l 

5c{b-l) 



Fig. 6 shows the monopole of the galaxy power spectrum, 
using the full expression Eq. (105) and Eq. (106) with 
Eq. (112). They generally agree very well, apart from 
the lowest-redshift case. 

The effective /nl given by Eq. (112) can serve as a use- 
ful tool to forecast for a given survey whether relativistic 
effects become important: if a survey achieves a fore- 
casted precision on the local /nl of order f^^, then rel- 
ativistic effects are relevant. Fig. 7 shows /^^ for differ- 
ently biased tracers as a function of redshift when Q = 0. 
Typical values are around 0.2, and generally lower than 
0.5. Note that we have assumed the universal mass func- 
tion approach to estimate /^f,; if for some reason the 
tracer number density evolves very rapidly with redshift 
so that be becomes large, f^ would increase correspond- 
ingly. In general however, we expect that number to re- 
main less than 1. /^^ rises sharply for z ^- 0, because 
of the 1/x factor in C [Eq. (97)]. Fig. 8 shows /°| for 
galaxies of bias 6 = 2 with different luminosity function 
with slope of (Q = —1, 0.5, 1, 2). This figure indicates 
that /^L varies with Q, but does not exceed 2 except 
for low redshifts where, again, f^ diverges due to l/x- 
This divergence however, is removed completely for dif- 
fuse backgrounds where Q — I (dashed line). 



CONCLUSIONS 



Future galaxy surveys will measure the clustering of 
galaxies on scales approaching the horizon. Thus, it is 
necessary to embed the observed galaxy density in a rel- 
ativistic context. This problem has received considerable 
attention recently. In this paper, we have derived the 
observed galaxy density contrast Sg in terms of the in- 
trinsic galaxy overdensity Sg and metric perturbations in 
the synchronous-comoving gauge (Sec. II). By transform- 
ing to conformal-Ncwtonian gauge, we reach agreement 
with the results of Refs. [11, 13]. On the other hand, we 
find some minor disagreement with the expression of [9] , 
which can be traced back to a sign issue. In App. C, we 
also show that our formula for Sg reproduces the correct 
analytic result for the following five test cases: 1) pure 
spatial gauge mode, 2) zcro-wavenumbcr gauge mode, 3) 
perturbation to the expansion history, 4) spatial curva- 
ture, 5) Bianchi I cosmology. These tests cover most of 
the parameter space of actual cosmologies in a simplified 
setting, and thus lend credibility to our result. 

One necessary, further ingredient for a description of 
galaxy clustering on large scales is a physical, gauge- 
invariant definition of galaxy bias. We present a straight- 
forward physically motivated definition in Sec. III. This 
bias relation is easily seen to reduce to the standard 
linear bias relation in synchronous coordinates, where 
constant-time hypersurfaces coincide with constant-age 
hypersurfaces. The bias relation Eq. (80) can be seen 
as a proper generalization of the peak-background split 
bias. Using this result, we arrive at a simple expres- 
sion for the observed galaxy density perturbations in 
synchronous-comoving gauge [Eq. (96)], described by the 
(gauge-invariant) bias parameter b, and the redshift evo- 
lution of the tracer population [through be, Eq. (54)]. 

The recent study of Baldauf et al. [26] has also reached 
the same conclusion by constructing locally flat space- 
time coordinates around a freely-falling observer, where 
long-wavelength modes locally act as curvature. In that 
coordinate system, the galaxy number density is modu- 
lated not only by the local curvature but also the time 
difference between global time and local time whose ef- 
fect exactly coincides with our be in Eq. (80). In contrast, 
Ref. [9] adopted a bias relation in the constant-observed- 
redshift gauge, which leads to considerably different pre- 
dictions for the large-scale galaxy power spectrum. How- 
ever, this does not seem to be a physical description of 
galaxy bias, since the age of the Universe is not constant 
on constant-observed-redshift slices. 

We believe these results, together with other recent 
work [11, 13, 26], allow us to unambiguously predict 
the two-point statistics of large-scale structure tracers 
on large scales. 
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FIG. 6: Monopole galaxy power spectrum [Eq. (108)] from the full expression Eq. (105), and using the effective /nl approxi- 
mation, Eq. (106). We use the same set of bias and redshift as in Fig. 5 (and Q = 0). The effective /nl approximation is valid 
within the horizon (fc > aH), marked by a vertical line, while it breaks down on large scales where the x~ term in Eq. (109) 
is important. The breakdown is more apparent when B is large (upper two panels). 





FIG. 7: Effective /nl [Eq. (112)] as a function of redshift for 
different galaxy biases when Q = 0. 



FIG. 8: Same as Fig. 7, but for galaxy bias (b = 2) with 
different magnification Q — —1, 0.5, 1 and 2. 
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Appendix A: Metric variables and gauge 
transformations 



Let us consider a general scalar coordinate transforma- 
tion (T, L) 



x" -^ x" = x" + (T(x"), d'L{x°')) 



(Al) 



While true scalar quantities are invariant under such a 
coordinate change, perturbations around the background 
do change because the background is time-dependent. 

For example, consider a scalar function C(t, x) whose 
background value only depends on time ({t) . This is the 
case for all scalar functions in the homogeneous universe. 
As a scalar, the function ( docs not change under the 
coordinate transformation in Eq. (Al), and it is only the 
background time r, thus ({''') ^ that is changed. From the 
relation 



C = C(T)+<5C(T,x)=C(f) + (5C(T,x), 



(A2) 



we calculate the scalar perturbation 5(^ in the trans- 
formed coordinate in terms of the variables in the old 
coordinate as 



SC = 5C + C(t) - C(t + T)=dC- 



dr ' 



(A3) 



up to linear order in perturbations. By appyling 
Eq. (A3), we find that the matter density contrast and 
redshift perturbation transform as 



-'Tn — ^Jm 



dlnpr, 
dr 



■T = 5„, + 3aHT, 



5z = 6z + aHT, and 



5g — 6g ~ beUHT, 



(A4) 



where in the second line we have used that dz/dr — ~aH. 
Hence, the quantity 6m, ~ 3(52: remains invariant under 
gauge transformations. The scalar component of the pe- 
culiar velocity, diV = adx^ /dt, transforms as 



V = v + L'. 



(A5) 



We can calculate the gauge transformation of higher rank 
vectors/tensors in a similar way, by applying the appro- 
priate transformation law of the object. 



Consider a general FRW metric with scalar perturba- 
tions defined through 

ds^ ^ -a^{l + 2A)dT^ - 2a^B^,dTdx' 

+ a^ [(1 + 2D)d^j + 2E,j] dx'dx^ (A6) 



where 



E,, = ( d,d, - -<5,,V2 1 E. 



(A7) 



Then the metric perturbations transform as follows: 

A = A- aHT - T', 
B =B + L' -T, 

D - \V^E = D- \w^E - aHT, and 

E = E- L. (A8) 

Similarly, when the perturbed metric is defined through 

ds'^ = - a2(l + 2a)dT^ - 2a^P^,dTdx' 

+ a^ [(1 + 2ip)S,j + 27,ij] dx'dx^, (A9) 

perturbations transform as follows: 

d = a — aHT — T , 

$^(3 + L' -T, 

if = ip — aHT, and 

7=7-L. (AlO) 

Note that the spatial metric components in Eq. (A6) and 
Eq. (A9) are related through 



a = A and ip ^ D - -V^E. (All) 

o 



2. From synchronous-comoving to 
conformal-Newtonian gauge 

The two most commonly considered gauges in cosmol- 
ogy arc the conformal-Newtonian gauge, which is defined 
through 

B = E = ^ /3 = 7 = 0, (A12) 

using the convention Eq. (A6) and Eq. (A9), respectively, 
and the synchronous-comoving gauge, defined through 

A = B = v = ^ a = f3 = v = 0, (A13) 

in the same conventions. Here we explicitly give the 
transformation between these two gauges, using the gen- 
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eral expressions of the previous section. We find that 

= a - aHT - T' , 
= L' -T, 



D- -V^E = ip-aHT, 
o 

E = — L, and 
= v + L'. 



Solving these equations leads to the transformation from 
conformal-Newtonian to synchronous-comoving gauge, 

D = >f+ -V^ / drv + aHv and E = drv, (A15) 

and the corresponding inverse transformation, 
a= - aHE' - E" , 



ip = D- -V^E-aHE' and 
^ 3 

v^E'. 



(A16) 



Note that in Eq. (A15), we have the freedom to add an 
integration constant, which reflects the fact that D and E 
in synchronous gauge contain a spatial gauge mode. Such 
a gauge mode can be removed by introducing (j) = D — 
y'^E/3 and/or taking time derivatives of E. For example, 
Eq. (A16) contains only and E' , because the conformal- 
Newtonian gauge has no residual gauge freedom. 



3. Metric variables in synchronous gauge 

In this appendix, we show the relation among three 
common parametrizations of synchronous-comoving 
gauge, and how to relate all metric perturbations to the 
matter density perturbation in the adiabatic case. 

The convention used in this paper is 

%,(t, x) = a^{T) [2D{t, x.)d,j + 2E,,{t, x)] , (A17) 

where Eij is traceless and given in Eq. (A7). In Yoo [10], 
the spatial metric perturbation is defined as 

%,(t, x) = a^T) [2,/)(t, x)J,, + 29,a,7(T, x)] . (A18) 

Comparing Eqs. (A17,A18), we find the relations 

(/)(t,x) = D{t,x) - -V^£;(t,x) and 7(t,x) = E{t,x). 

(A19) 
In Fourier space, the metric in Eq. (A17) becomes 



(55,,(T,k) = a^T) 2D{T,k)S^ 



1 



-2[Kk,-^e6^j]E{T,k) ,(A20) 



which can be compared with the spatial metric pertur- 
bations defined in Ma & Bertschinger [31]: 



Sgij{k,T) = a'^ir) 



+ 6 






p - g-^'J ]vik,T] 



(A21) 



(A14) Comparing Eq. (A20) and Eq. (A21) leads to 



Dir, k) = ^^(^ and E{r, k) = -^("' ^^ + ^^("' ^^ 



6 



2fc2 



(A22) 



hJT, k) + 6??(t, k) 
2fc2 



and combining above results yields 

(/)(t, k) = — 77(t, k) and 7(r, k) = - 

(A23) 

We now relate the metric variables to the matter den- 
sity contrast in synchronous gauge. These relations as- 
sume General Relativity, Gaussianity, and adiabaticity, 
and are only used for the equatiosn and numerical re- 
sults in Sec. IV. First, let us consider the time derivative 
of E. The continuity equation in synchronous comoving 
gauge is given by [31] 



SUr,k) = -^h'iT,k), 



(A24) 



and from the Einstein equations we have [e.g. Eq. (22) 
of [31]] 

k'^v'ir, k) = AinGa'^k^T^ (xv^O. (A25) 

Therefore, we calculate the time derivative of E as 

(A26) 
where / = dlnD/dlna. From here, we can also calculate 
the second derivative of E as 



S"(r,k) 



1 d[aHf6^ir,k)] 



k^ 



^2 rr2 



dr 



<5,„(T,k), (A27) 



where the second equality comes from the time evolution 
of the linear density contrast (continuity and Euler equa- 
tions). Finally, from the Einstein equation [Eq. (21a) in 
Rcf. [31]], we calculate (^(r, k) = — 7y(T, k) as 

<^(r,k) = - ^ [aHh'{T,k) + SnGa^ST^] 



{a'H'f + A^Gp^a') 



2,(5m(r,k) 



a2^2 J ^ ^^^ 



fc2 
Srn{T,k) 



(A28) 



Note that while D contains a spatial gauge mode, cj) = 
D — V2_B/3 does not contain any. Also, from Eq. (A25), 
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it is obvious that (/) is constant in time. The physical 
interpretation is that for a plane wave perturbation, two 
neighboring test particles separated by an infinitesimal 
distance perpendicular to k have a separation that is pro- 
portional to the background a{t); only the component of 
separation parallel to k is perturbed. 



Appendix B: Derivation of Eq. (96) 



The second, integral term is clearly equal to the ISW 
contribution, and neglecting it we obtain Sz = [d\\E' + 
E"]l, which in Fourier space becomes 



5z = ik^iE' + E" 



(B6) 



Here, we have ignored the contribution from the origin, 
and jjL denotes the directional cosine between the wave 
vector and the line of sight direction. 



This section outlines the derivation of Eq. (96) from 
the expression for 5g [Eq. (55)]. We begin by deriving 
a compact expression for (5z, before moving on to the 
convergence k and the derivation of Eq. (96). 



1. Integrated Sachs- Wolfe term in synchronous 
comoving gauge 

Here we identify the ISW term in synchronous- 
comoving gauge and derive a simplified expression for 5z. 
In conformal Newtonian gauge, the ISW term is given by 



iJzisw 



/ dx{^ + -^)', (Bl) 

Jo 



where $ and $ are the Bardeen potentials and prime 
denotes a derivative with respect to the conformal time 
T. In synchronous comoving gauge, they are 






aHE'-E" and 
D + -V^E + aHE', 



and the ISW term becomes 



Szisw = f dx (d' - iv^S' + E" 



(B2) 



(B3) 



Using Eq. (30) and the definition of £', we find that the 
redshift perturbation in the same gauge is given by 



5z 



dx ( D-\v'^E + dlE 



(B4) 



Now, we use that d^ = d\\ — 9r, and rewrite the red- 
shift perturbation by successively applying integration 
by parts: 



6z 



dx 



D' -]-V^E' + {d^+dr)d\\E' 



dx 



D'--W^E' + (d^ + dr)E" 



[d\\E' + eX + £ dxl^D' -^V'E' + E 



(B5) 



2. Convergence in synchronous comoving gauge 

In this section, we calculate the coordinate convergence 
in the synchronous-comoving gauge defined in Eq. (48). 
By applying the perpendicular derivative d±i to the per- 
pendicular directional displacement Ax'\_ in Eq. (39), we 
find 



1 /"X ~ 

^= -o^i / dx{x-x)-{D + E\\) 

2 Jo X 

+ d^, dx- (E]fi^ - E\\k 

- \^^^ [x[E^Ao)fl^ - Ejkiopn^n'^]] 



(B7) 



where we have labelled the three terms and used "(o)" to 
denote metric shear evaluated at the observer. In the first 
term, we have pulled out the perpendicular derivative 
inside the integral over the unperturbed geodesic. Hence 
the additional factor of x/x in the integrand. 
We begin with the ky^> term: 



x(i) 



ivi j^ dxix -x)^{d- I'^'E + d^E 



K V 



^-' ^''dx{x-x)-{dlE~E") 
X " 



(B8) 



This allows us to relate k to k with the introduction of 
some new terms. We can combine these terms with terms 
in kP'' if we simplify the latter: using the commutation 
relations Eqs. (9)-(13), we have 



n^E'^-f^Ew 



n' { d,d' ~ U'jvA E - n' (df^ - ^\7^ ] E 



= (n^djd' - ed^] E = (d\\d' - d\\edA E 



d«d]E^d]_d«E--d]_E, 

X 



(B9) 



and so 



''"^^'^£HP''-P)^ '"'»' 
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(Note again the additional factor of xlx that arises since 
when we move d\ outside the integral it acts at radius x 
rather than x.) Combining these gives 



/j(l)+^(2) ^ ^ 



1 . 



LV X 



-2 ~2 

1 /■* 
^+2^1 J dx 






»);|l2^.^- 



dE 
dx 



-9 -9 

V V 



(Bll) 

where in the second line we have used ' = d\\ — d/dx- 
The duE terms form a total derivative, which may be 
separately evaluated: 

d^E 






X' - 

X 

X 



dx^ 



-vi 



x! 

X 



X 5||£; 



2^E 
X 

(B12) 



The quantity in braces vanishes at x, whereas at e — >■ it 
blows up. For this reason, we will evaluate this expression 
only for e > and then take the limit after all divergences 
are cancelled. If we Taylor expand d\\E to order e, we find 

d\iE{e) = fi'E,{o) + fi'k^eE,j{o) - h'E\{o)e + 0{t^). 

(B13) 
The quantity in braces is then (keeping terms that are 
nonvanishing as e — )■ 0+) 

- f- Y + x) E^fi' + x^fefi^E^,,{o) - x^ffE^io). 

(B14) 
The V^^ operator pulls down a factor of for monopoles, 
— 2/x^ for dipoles, and — 6/x^ for quadrupoles. It follows 
that 



:-(i) 



^(2) = ^_ 



dx 



X^ \d^E x^ 

— -xhrT-2^^ 

X ^ ax X 



X/ 



e 



-2n'£;',(o). 



The remaining integral is also a total derivative, although 
in this case two integrations by parts are necessary. Both 
terms can be integrated by parts to obtain a single deriva- 
tive dE/dx', the integrals cancel leaving only the bound- 
ary terms: 



AiW+/i(2) ^ ^. 



iv^ 



x! 

X 



\dE ^ 



X 



+ [-- + i] E^,(o)n' - QE^j{o)n'n^ 



e 



-2n'E'^{o). 



(B16) 



This simplifies to 



«(i)+^(2) 



-+2Vi 



X' ~\dE,. x\.. 



1. 



:Vi£; 



e X 



E4o)h' 



-QE,j{o)fefi^ +2fi'E'/o) 



(B17) 



where V5_ is evaluated at radius x- 

Further simplification requires the limiting forms of the 
terms in Eq. (B17). A lowest-order expansion gives the 
quantity in brackets as 



x' 



X [n'E,{o) - E'{o) + ek'fi^E,,{o) 



-2en'E[^{o) + eE"{o)] 

-\[E{o)+€n'E,{o)^eE'{o) 

+ ^e^fi^h^E^,,io) ~ e^fi^E^io) + ^e^E"io)]. (B18) 

Again using that V^^ operator pulls down a factor of 
for monopoles, — 2/x^ for dipoles, and — 6/x^ for 
quadrupoles, we find that Eq. (B17) simplifies to 



Ai(i)+«:(2) 



K+-ViE+-Ei{o)h' 
2 X 



--E,,io)fffi^~eEr{o). (B19) 

The third term, k^^', is independent of x since it is a 
derivative of a quantity that depends linearly on x; hence 
we may evaluate it on the unit sphere x = 1- We find 



^(3) 



--E,j{o){6ij - n'fi^) + -Ejk{o)[2fi^fi'' 






f)*Ti'^~lf)*Ti-Jl 



3 



-E,j{o)ev. 



(B20) 



Combining with Eq. (B19) gives 



k^ K+ -W]^E + -Ei{o)il' - h'E'Ao). (B21) 

2 X ' 



(B15) This is Eq. (51). 



3. From Eq. (53) to Eq. (55) and Eq. (96) 

Let us start from Eq. (53), 

Sg(Si.) =dg + be5z + 2D-E\\ 

Aa:|| / l + SdH(S)\ , 

+ 2—r^ -2k- { I ^ 1 Sz 

X \ H dz 



ig(^'+^i) 



(B22) 
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where 



Aa;||=-/ dx{D + E^) - ^:^5z 



His) 



(B23) 



and we use Eq. (B21) for k. Also, in the previous sec- 
tion, wc found that the rcdshift perturbation, ignoring 
the ISW term, is given by 



Sz^diiE' + E" 



(B24) 



As shown in App. A2, _D and E contain spatial gauge 
modes, while <j> = D — V^i?/3 and E' remove such gauge 
modes. We now collect the four terms in Eq. (B22) which 
contain gauge modes: 



2D 



AxH 

En +2 ^ 

X 



2k 



20 + ^S/^E - (^df iv^) i? - I £ dx{D + i^ii 

- -^-^Sz -2k- \I\E - -fi'E.io) + 2fiE'^{o) 
xH(z) X ' 



20- 



l[dnE-dnEio)]-i f dx{^ + d?E) 



X 
21 + z ^ 



xH{~z) 



—5z-2n + 2d»E'{o). 



(B25) 



The third term can be further simplified by double inte- 
gration by parts to yield 



dx{4> + dlE) 



[ dxi^ + E")+E' + d\\E 
Jo 

-E'{o)-d\iE{o), (B26) 

Then we find for the observed galaxy density contrast 

Sg{x) = Sg + beSz - 

1 + 5 dH(S) 2 1 



1- 



H{~z) II 
5 dH{2) 



H 



dz 



xH{z) 



Sz + 2(t) 



2K + 2duE'{o). 



XJo 



dx {<t> + E") 
(B27) 



This is Eq. (55). Note that the magnification bias contri- 
bution Q5M does not contain gauge modes [Eq. (72)]. 
To proceed to Eq. (96), we drop several terms: 

• The observer terms (which contribute only to the 
monopole and dipole). 

• The 0' term, since the Oi component of the Einstein 
equation ensures that in a ACDM universe 0' = 

[see Eq. (A25) and note that ry = —0]. 

• The terms -I j^ dxifp + E") in Eq. (55) and 
Eq. (72), corresponding to the time delay which 
is very small [14, 32]. 



• The terms involving the convergence k, which is 
generally not small but is a projected quantity, and 
in the flat-sky limit contributes only to transverse 
modes. 

Substituting in Eq. (B24) for 5z and using Eq. (72) 
then yields Eq. (96). Conversion to Fourier space in the 
flat-sky limit - i.e. where we make the replacement d\\ — > 
ifc/i - then gives Eq. (99). 



Appendix C: Test cases for the observed galaxy 
overdensity 

This appendix considers several analytical test cases 
that serve as a cross-check of Eq. (55). The first two 
cases are pure gauge modes, with the expected result 
that 5g does not receive any contributions from such per- 
turbations. We then consider a perturbed expansion his- 
tory, spatial curvature, and a Bianchi type I cosmology 
with anisotropic expansion. Finally we consider a model 
with a time-dependent linear gradient in E^ which has no 
metric perturbation but leads to a delicate cancellation 
of terms in the galaxy density. In all cases, the linearized 
version of the exact result can be derived straightfor- 
wardly, and we show that it agrees with the prediction 
of Eq. (55) in all cases. We further evaluate the magni- 
fication M using Eq. (72) and show that it matches the 
expected results. 



1. Pure spatial gauge mode 

The residual gauge freedom of the synchronous gauge 
allows us to reparameterize the spatial coordinates ac- 
cording to x* —> x* + ^' , where ^* depends only on the 
spatial coordinates and not x". This leads to a spatial 
metric perturbation hij = ~a^(C* j +^"',i)- Since here we 
consider scalar perturbations only, ^' should be derived 
from a potential ^' = ^^i. If we start from an unperturbed 
Universe, the resulting metric perturbation is 



D = — V^C and E = -£,. 



(CI) 



Eq. (CI) corresponds to a pure gauge mode. 

For this mode, we proceed to evaluate Sz, (p, and k. 
Since D and E are time-independent, it is trivially seen 
that 5z = 0, and Eq. (CI) immediately implies = 0. 
Finally, since D = \^^E and E" = we can also see 
that K — {). Then Eq. (55) reduces to 



4w 



(C2) 



which is the expected answer. That is, in this case we do 
not have any contributions to the galaxy density aside 
from the intrinsic contribution. 

Wc may also evaluate the magnification using Eq. (72). 
With E' = Q everywhere and 5z = (f) — n = Q, it is 
trivially seen that SM = 0. 
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2. Zero-wavenumber gauge mode 

There is another spatial gauge mode that does not fall 
into the rubric of Eq. (CI): the zero-wavenumber mode 
given by 



Z? = 5 + gix'' and E = 0, 



(C3) 



where S is a constant scalar and g is a constant vector 
field. This is generated by the gauge perturbation 



Again we trivially have dz = 0, but this time 
Xg ■ n. Also we have 



-2 



(C4) 



Vl7^ = Vi(xg-n) = — g-fi, 



X 



(C5) 



since g • n is a dipole {£ = 1) and for a pure multipole of 
order £ the operator V^ yields a factor of —£{£ + l)/x^. 
Thus we find 

K = --/ dxix-x)- — g-n=-xg-n. (C6) 
2 Jo XX 2 

The galaxy density perturbation obtained via Eq. (55) 
has only four nontrivial terms: 

Sg = Sg+2{^+Xg-n}-- / (^+xg-n)dx-Xg-n. (C7) 
X Jo 

Here the second term comes from the +20 term in 
Eq. (55), the third term is the line of sight integral of 
(j> + E" ~ (j>, and the last term comes from the —2k. It is 
easily seen that these three terms cancel, leaving 6g ~ 6g, 
which is the expected answer. 

Unlike the previous case, here the magnification con- 
tains nontrivial terms: substituting the nonzero values of 
(j) and K into Eq. (72), we find 

SM = -2(E + xg-n) + xg-n+- / {E + xg-n) dx = 0, 

xJo 

(C8) 

as expected. 



3. Perturbation to the expansion history 

A less trivial type of perturbation is one in which we 
alter the cosmic expansion rate. This can be done by 
setting 



D = D{t) and ^ = 0. 



(C9) 



The Universe so described is still an FRW model since 
it is homogeneous and isotropic. (It may no longer be a 
solution to the Friedmann equation with only matter-|-A, 
however this does not concern us since we are testing 



an equation derived only using kinematics.) However, 
it has a "true" scale factor atruc that is related to the 
unperturbed scale factor via 

atruc(r) - a(r)[l + D{t) - D{to)], (CIO) 

where we fix atruc to be unity today. The true time co- 
ordinate (proper time in the case of FRW) remains equal 
to the coordinate time. 



The true conformal time is then 
rf^tnio [ adr 

Ttr 



itruo = t= / a{T)dT. (Cll) 

hen 

= [[l-DiT)+D{To)]dT. 



^tiuc J ^tru 

(C12) 
(The integration constant is chosen to set rtiuc = at the 
Big Bang, but we do not need to make use of this fact.) 
Integrating gives 

/To 
Z?(Ti)dTi. (C13) 

We care in particular about the behavior as a function 
of the observed redshift z, which is related to atruc = 
(1 -I- z)~^. It follows from Eq. (CIO) that the comoving 
distance relation is now 



air) = at,-ne{T)[l ~ Dir) + D{to)] 



(C14) 



and so we may write the perturbation to the conformal 
time. 



T (atruc) = Tbg(atruc) + 



D{to) - D{t) 



(C15) 



atruc^(atruc) 

Here Tbg is the background conformal time-scale factor 
relation. Finally using Eq. (C13) yields 

Ttruc.O — '''true = [1 + -D(to)] [^bg (1) - ^bg (fltruc)] 

D{t^) - D{t) 



dtrueH {atr 



/To 
D{Tl)dTl. 

(C16) 



This is the true comoving radial distance Xtrue to redshift 
a^ig — 1. That is, 

Xtruc = [1 + -D(To)]Xbg(atruc) 

DJTo) - D{t) rro 

O-tiucH [atruc) 

The true Hubble rate at this time is 
_ dlnatruc 

-Htruc — 



D{Ti)dTi. (C17) 



dt 

H{r)- 



D'{t) 



^ dH. . , , ,. D'{t) 

^bg(atrucj + -3— F(atruc) - '''bg(atruc) H 

dr a 



^bg(atruc) 



dH D{To)-D{T) ^ D'{t) 
dr atruc-ff(atruc) a 

dH,^, . „, ,, D'{t) 



Hbgiatruc) + a—[D{To) - D{t)] + 



(C18) 
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We expect the perturbation in the observed galaxy 
density to have several parts: there is a perturbation 
in the physical galaxy density, a part associated with 
the different epoch in cosmic history at which the galaxy 
density is measured (different t; one wants the different 
physical density here so we include both the be evolution 
and the —3 associated with the dilution of comoving vol- 
ume), and a part associated with the different physical 
volume. Specifically: 



dg^Sg + {he - i)[D{To) - D{t)] + In 



dVtmc I datmc dil 
dV/dadn 



(C19) 

The comoving volume effect is computable from 
Eq. (C17). We see that 

dVtruc ^trucXtruo 



"fltruc dll Hfruc 

.2 



atruc[Xbg(atruc)]^ ( ^ r,jj^ 
HhJatruc) ^ 



to) 



2j;°D{n)dT, 2[DiTo) - D{t)] 



X(otrue) atruoi?(atiuo)x(atrue) 

r a dH ,^, , , ,, D'{t) i 



} 



(C20) 



This leads to 



5g - 5g + {be-?,)[D{T^)-D{T)]+2D{T,,) 



dlnH 



d\na aHx 



[D{to) - Dir)] 



2pD{T,)dT^ D'jr) 
X(atruo) aH 



(C21) 



The expected "magnification" 5M is twice the perturba- 
tion to the angular diameter distance at fixed observed 
redshift; in a flat universe this is equivalent to the per- 
turbation to Xtruo- Using Eq. (C17), 



A simple comparison shows this to be equivalent to 
Eq. (C21). Similarly, evaluation of Eq. (72) gives 

5M = -2D{t) + - / D{Ti)dTi 
xJt 



aHx 
in agreement with Eq. (C22). 



[D{to) - D{t)], (C24) 



4. Spatial curvature 

A fourth example of a perturbation we consider is 
spatial curvature. Under stereographic projection, a 3- 
sphere of curvature K (radius of curvature K~^'^) can 
be written with 3-metric 



1 

4' 



dsl= il + 4^|xp 1 dx'dx\ (C25) 



or to first order in K, 



D = — Alxp and E = Q. 
4 



(C26) 



The expected result is that in the perturbed universe, 
the radial comoving distance-redshift relation remains 
the same. However, there is a change in the volume ele- 
ment associated with the change in the comoving angular 
diameter distance. 



rfKow _ sin|-x 

■02 



= 1 - Ikx^ + 0{K^), (C27) 



dVoiA X' ' 3 

where siui^ is the sinelike function: 

X K = Q 

siTiK X = { K"^/^ sm{K^/^x) K >0 

(-A:)-i/2 sinh[(-A')i/2;^] K <0. 



Thus we expect to obtain 



(C28) 



5M 



, Xtrue Xbg 



Xbg 



-2D{to) + 



2[D{to) - D{t)] 



aHxhg 

/ DiT,)dT^. 

Xbg Jt 



(C22) 



In comparison, if we use Eq. (55) to find 6g, then we 
find that the perturbation in Eq. (C9) yields <f> = -^'(r), 
K = 0, and 6z — D{to) — D{t). Therefore, 



Sn = S„ 



6.-1- 



d\nH 



^2D{t) 



d\na aHx 
2pD{T,)dT, D'{t) 
X 



[D{t,)-D{t)] 
(C23) 



aH 



h^^g-^Kx"- 



(C29) 



The magnification 5Jv[ is —2 times the fractional per- 
turbation to the angular diameter distance coming from 
spatial curvature, which is 



5M 



,sin^^|-x ^ 1^.^.^ 



(C30) 



If we instead use Eq. (55), we find that = — ^A'x^, 
K = (since D is a pure monopolc, V\D = even 
though V^Z) ^ 0), and 5z = 0. Then 



l^-~2 2 fV l^. 



Evaluation of the integral trivially recovers Eq. (C29). 
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We can also compute the magnification from Eq. (72); 
we get 

SM = ^Kx^ + jp ("i^^') '^^ = 5^^'' (^^2) 
in agreement witli Eq. (C30). 

5. Bianchi I cosmology 

The previous test cases have not tested the terms in- 
volving £". One case that does is the Bianchi I cosmol- 
ogy, in which the three spatial axes (usually taken to be 
the coordinate axes) have different scale factors but the 
universe is still homogeneous. We will focus here on the 
case where the observer looks in the x'^-direction and the 
metric perturbations are 



£'(x, t) = — S3(t) and 

tj[yi,T) = 



,(C33) 



with si(t) + S2(t) + s^{t) = and Si(To) = 0. This is 
equivalent to a case where the expansion along the 3-axis 
is unperturbed, but the other two axes are perturbed: the 
scale factors are 



ai{t) = a{t)[l 
a^it) = a{t)[l 
a3.{t) = a{t). 



Si(t)-S3(t)], 

S2{t) - sz{t)], and 



(C34) 



Note that we have already considered in Eq. (C9) the case 
where the global isotropic expansion D{t) is perturbed, 
so no new independent tests of our results are possible 
by using a different function for D in Eq. (C33). Also 
we have considered in Eq. (CI) the case where ii^ is a 
time-independent function with zero Laplacian, so there 
is no independent test of our result that can be obtained 
by allowing Si{ta) ^ 0. 

It is straightforward to determine the expected change 
in observed galaxy density for this model. The metric 
components {30075037533} are not perturbed, so for the 
n = (0, 0, 1) direction the past light cone is unperturbed 
and a given redshift z corresponds to the usual distance 
Xbg('Z) and Tbg(2:). The only nontrivial effect is in the 
transverse dimensions and in the volume element - the 
angular diameter distance is modified by the perturbed 
expansion rates in the 1 and 2 directions. We may de- 
termine the true angular diameter distance along the 1 
axis by considering a ray projected backward from the 
observer with a physical angular separation <; from the 
3-axis, i.e. in direction n = (<^, 0, 1). (We work to order <^ 
so that sine; = <;^ and cos<;^ = !•) Then the 4-momentum 
of such a ray with unit energy is 



V^. = {-l,-^.Q,-l) 



(C35) 



Since the metric coefficients do not depend on spatial 
position in this model, the spatial covariant components 



Pi of the momentum are conserved. Then we find that 
the spatial position is given by 



dx^/d\ 



dri 



' [ai(n)]-V 



'TO [a(n)] % 
= -^ / [l-2si(Ti) + 2s3(Ti)]dri. (C36) 

J Tq 

The physical angular diameter distance is the physical 
transverse distance divided by the angle subtended, i.e. 

a\[T)x^ l^,. That is, 

-DA.phys.l = -a(T)[l + si(t) - S3(t)] 

X / [l-2si(Ti) + 2s3(Ti)]dri 

J To 

= a(T)|l + Sl(T)-S3(T) 

^,r[^3(n)-.r(n)]dny ^^^^^ 

To - T J 

The observed galaxy overdensity then deviates from the 
true galaxy overdensity only by the transverse area ele- 
ment (since the time of observation and the longitudinal 
distance-redshift relation are unaffected). That is. 



'^9 = ^9 + 1^ 



£'A,phys,l£'A,phys,2 

n2 ■ 

A,unpcrt 



(C38) 



Using that X]i=i ^i(''') = 0, we may simplify this to 



4 = '^9 ^ 3s3(r) +6 



\l° S-i(TY)dTY 
Tq-T 



(C39) 



The magnification should also be given by the change in 
the transverse area element: 



bM = 3s3(r) -6 



To -T 



(C40) 



We now wish to compare our result to Eq. (55). By 
construction we have on our chosen sightlinc _D -I- £'|| = 
0, so 6z = 0; and since S/'^E = we have cj) = —S3. 
The convergence k is more complicated: D is spatially 
constant, V^-E = 0, and E is, a, pure quadrupole {i = 2) 
and hence V^^ pulls down a factor of —£{(, 4- l)/x^ = 
— 6/x^7 so wc find 



V2 [D--y^E-E"\ = -%- 

= -34', 



6(4'xV2) 



(C41) 



where the second equality is valid only on the 3-axis line 
of sight. Consequently the convergence is 



K= - I dx{x-x)-S3ix)- 



(C42) 
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Finally we have d'nE' = Sg, and the observer terms E'{o) 
and 9||£"(o) both vanish. Plugging these results into 
Eq. (55) gives 



1 



xJo 



■4(x)-2s3(x)-^ 



X 



2X^4(x) 



1 



-^3(X)+2XM(X) 



-3 

/o 

l + z 

H 



dx{x~x)-4ix) 

A 

-4(x)]- 



(C43) 



The terms containing s'^{x)/ H cancel and the two inte- 
grals can be combined, yielding 



5,-2,S3(x)-x4(x) 



dx 



|.3(x)+(2^-3x)4'(x) 



.(C44) 



This does not quite resemble Eq. (C39), but we can cast 
it in a similar form by applying repeated integration by 
parts to the second derivative term. For a general func- 
tion /, 



/(x)4'(x)rfx = -/(x)4(x) + /(o)4(o) 



+ /'(X).S3(X)-/'(0)S3(0) 

+ r f"{x)s,{x)dx, (C45) 



where the unusual signs result from the fact that ' denotes 
a derivative with respect to r instead of x (the relation is 
simply a minus sign). Then, recalling that sz{x = 0) = 0, 
Eq. (C44) simplifies to 



5g^5g- 3s3(x) + dx -S3(x)- 



6 

X 



(C46) 



Inspection shows that this is equivalent to Eq. (C39) via 
a change of variable, t = tq — x. 

For the magnification, Eq. (72) predicts 

SM = 2s3(x)+x4(x) + 3 r(x-x)^4'(x)rfx 

Jo 



2 fx 
xJo 



-ssix) + 2X'4'(X) 



X 

dx; 



(C47) 



repeated integration by parts again reduces this to a form 
equivalent to Eq. (C40). 



6. Potential-only mode 

The only terms left in Eq. (55) that we have not tested 
are the observer terms, E'{o) and d^^E'{o). These can be 
tested using a potential-only mode 

D{x,t)=0 and E{yi,T) ^ T{t) + w,{t) x\ (C48) 



This mode has no metric perturbation, D — Ei- 



0, 



and so we expect to get 5g ~ Sg. However it does have 
nonzero observer terms. 

Trivial evaluation shows that for the "perturbation" 
Eq. (C48), we have Sz — (f) — 0. However the convergence 
K is not zero despite the vanishing metric perturbations! 
Instead we have 



1 



'dx(x-x)^Vi(T" + <n'x) 
X 



dx[i-4] w['fe 



(C49) 



where in the second line the action of V^ is to eliminate 
the monopole and extract a factor of — 2%"^ from the 
dipole. Now since Wi is a function only of r we have 
w'l = (f'Wi/dx^ . Double integration by parts then gives 

(C50) 

X 

Similarly we find 



-I II N .jU'j(To)-Wi(ro-x) 

X 



dxE" 



dx (T" + h'w'/x) 

T'{ro)-T'{ro-x)-xn'w',{To-x) 
+n'-Wi{To) -n^Wi{To- x)- (C51) 



Finally the remaining terms are 

E'-E'{o)^T'{To-x)+xh'w',{To-x)-T'{To). (C52) 

and 

d\iE'{o)^n'w'^{To). (C53) 

Assembling these pieces of Eq. (55) then leads to a mass 
cancellation that recovers Sg = 6g, as expected. 

The magnification equation, Eq. (72), has the same 
nonzero pieces and a similar mass cancellation occurs, 
leaving the correct result SAi = 0. 



Appendix D: Connection with results in the 
literature 



1. Yoo et al. 



In this section, we compare our result, Eq. (53), 

~6g{i) = 6g - (l + S)- 



dhi{a^ng 



1 



dz 
l + zdH 
H H 



Sz 



Sz + 2D- E\\ 



-f 2 ^ -2k 

X 



X 

(Dl) 
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to Eq. (36) of Yoo et al. [9], restricted to synchronous- 
comoving gauge: 



b{6m- iSz) + 2D 



Eiju'h^ 



1 + z — <5z-2— — 
oz Hr 



Sz — Sz 



- 5pSDl - 2k 



1 + SdH 
H "d7 



■6z- 



St 
2—. (D2) 
r 



We can convert their result to our notation by noting 
that r ~ X &nd 



Sr 



1 + ^ 
HCz)' 



(D3) 



Note that the Sr defined in Eq. (16) of Yoo et al. [9] 
has a different sign of Eu compared to ours, as discussed 
in Sec. II. Further, d6z/dz = l/Hd6z/dx, and using 
Eq. (30), Eq. (D2) becomes 



Sohs =b {5„ 



1 

Ax|| 
X 



35z) + 2D 
1 + zdH 



H dz 
- ^p5Dl 



-E\\ 
5z- 

-2k. 



H{~z) ^ 



E[: 



(D4) 



Comparing the two expressions, we find two differences: 
first, our result involves the time-dependence of the num- 
ber density of tracers, d\afig/dz, while this quantity 
does not enter Eq. (D4) since the bias is defined in the 
uniform- redshift gauge (see Sec. III). The second differ- 
ence is the sign of the i?|[— term, which goes back to the 
difference in sign in 5r [our Eq. (38), and Eq. (16) in Yoo 
et al. [9]]. This was discussed in Sec. II. 



2. Challinor & Lewis and Bonvin &i Durrer 

We now transform our Eq. (55) into variables in con- 
formal Newtonian gauge in order to compare our re- 
sults with [11] and [13]. By using Eq. (A4), Eq. (A5), 
and Eq. (AlO), we first find the transformation law 
from synchronous-comoving gauge to conformal Newto- 
nian gauge as 

* = V^^^^ =a= -E" - -HE', 
$ = (/,(CL) ^_(^^ -(f> + HE\ 
v=E', 
Sz=6z + nE', 
Sm. ^5,n + 3nE', and 
Sg ^ 5g - be-HE' , (D5) 

where the " symbol denotes the quantities in the con- 
formal Newtonian gauge [Eq. (A9)], and the superscript 
(CL) denotes the metric perturbation variable defined in 
[11]. Bonvin and Durrer [13] use the Bardeen potentials 
^, $ as defined in the above equations. Note that we use 



H = aH in order to facilitate the comparison. We now 
transform the terms in Eq. (55) as follows: 

Sg + beSz =5g + beSz 

H dz xH{z)J 1-H2 xn r ' 



20 = 


^2{ip + 'Hv) 


X 


2, 
X 




-E" - 


= -{a- 


-^). 



(D6) 

By using Einstein's equations, i.e. (()' ~ Q [Eq. (A25)], 
and the definition of k 

^ = - o / dxix-x) ?Vi {a-^) = K, (D7) 

2 Jo X 

which is the same as the convergence in conformal New- 
tonian gauge {k defined in App. B2), Eq. (55) becomes 

~8g{k) ^5g + bjz - ^dlv - f ^ + A j [s, _ uv] 

2 2 f^ 
+ 2[f> + 'Hv]- -i!- - / dx{a-Lp)-2k 
X xJo 

= (5„ + bJz d'^v - I ^ H I 5z 

^ -H II yn^ xH) 

(i 2 /■* 

+ 2(p + 3Hv + a — —- — z dxia — (p) — 2k, 
'H xJo 

(D8) 

where in the second equality, we used the identity 

^v = a + Hv-^, (D9) 

and we drop quantities at the observer's position as is 
done in [11]. 

The redshift perturbation Sz in conformal Newtonian 
gauge can be calculated from Eq. (B5) and the gauge 
transformation Eq. (D5) as 



diiv — a — / dx{a — (p) 



(DIO) 



Rewriting the equation in terms of the variables in Challi- 
nor and Lewis [11] (a ~ tp, (p = —0, and (9||?) = v • fi), we 
find the redshift perturbation as 

<5z = V n - ,/.(CL) _ rdx^^^^'^'^+c^'^'''''^)' ■ (Dll) 



A similar change of variables can be made for the ob- 
served galaxy density contrast in Eq. (D8), which yields 



5g(i)=A„(n,z) + 3Hi; 



(D12) 
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where A„(n, z) is defined in Eq. (30) of Challinor and 
Lewis [11]. The additional term of 3Hv is due to the 
fact that the overdensity A„ is defined with respect to 
the physical, rather than comoving galaxy density. In 
performing the gauge transformation i — >■ i = i + T, we 
thus obtain an addition term from Ino'' = Ina'^ + 3aHT = 
lna3 + SHv. 

To compare the magnification terms, we transform 
Eq. (72) to the conformal Newtonian gauge, yielding 

2 2 f^ 

5M=- 2{ip + 'Hv) + -v + 2k- - dxia- ^) 
X xJo 

-2(i-^)is.-m) 

^2k + 20(CL) -I r dx f 0(^^) + ^(^^' 
X Jo ^ 

-^ -' (D13) 



'Km-''''^ 



which has to be compared to all terms oc 5s in Eq. (37) 
of [11]: 



K — r / dx{4> + i)) + 4) 
xJ 



-i' - / dx(0 +ip) + n-v 



(D14) 



As 5s is the same as 2 Q in our notation, our formula for 
6Ai is also in agreement with [11]. 

We can also write Eq. (D8) in terms of the variables 
in Bonvin and Durrer [13], 



Sg + beSz =Ds, 
a =*, 



9||W = — V • n, and 
djv =dr (V • n) , 



(D15) 



to arrive at the same relation Eq. (D12), after identifying 



rs -r 



K = - — / dr 



and A„ = A as defined in Ref. [131 



($ + *) 



(D16) 
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